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ABSTRACT 

In this article we identify and discuss various statistical and systematic effects influencing the 
astrometric accuracy achievable with MICADO, the near-infrared imaging camera proposed 
for the 42-metre European Extremely Large Telescope (E-ELT). These effects are instrumen- 
tal (e.g. geometric distortion), atmospheric (e.g. chromatic differential refraction), and astro- 
nomical (reference source selection). We find that there are several phenomena having impact 
on ~100yL(as scales, meaning they can be substantially larger than the theoretical statistical 
astrometric accuracy of an optical/NIR 42m-telescope. Depending on type, these effects need 
to be controlled via dedicated instrumental design properties or via dedicated calibration pro- 
cedures. We conclude that if this is done properly, astrometric accuracies of 40/ias or better 
- with 40/ias/yr in proper motions corresponding to x20 km/s at 100 kpc distance - can be 
achieved in one epoch of actual observations. 

Key words: Telescope — Astrometry — Instrumentation: high angular resolution — Tech- 
niques: high angular resolution 



1 INTRODUCTION 

The future optical/near-infrared European Extremely Large Tele- 
scope (E-ELT; see, e.g., Gilmozzi & Spyromilio ( 12008b '), which is 
designed with a 42-metre aperture, will offer a substantial improve- 
ment in angular resolution compared to existing facilities. At wave- 
lengths A = 2fim, diffraction-limited resolutions of © lOmas will 
be achieved. In terms of angular resolution in the near infrared, 
the E-ELT will outperform existing 8~10m-class telescopes like 
the VLT or Keck by factors of x4~5 and the future James Webb 
Space Telescope (JWST) by factors of »?. This increase in angu- 
lar resolution should translate into a corresponding improvement in 
astrometric accuracy. 

In order to exploit the E-ELT's resolution, a German-Dutch- 
Italian-French consortiurrQ proposed the Multi-AO Imaging Cam- 
era for Deep Observations (MICADO) in February 2008. As the 
spatial resolution of any ground-based observatory is initially lim- 
ited by the atmospheric seeing, MICADO will be equipped with a 
multi-conjugate adaptive optics (MCAO) system for achieving the 
diffraction limit of the 42m-telescope. This system uses three nat- 
ural and six laser guide stars for correcting the atmospheric turbu- 
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lence in a wide (> 2') field of view (Diolaiti et al. (2WSS). Images 
will be recorded by an array of 4x4 near-infrared (NIR) HAWAII- 
4RG detectors with 4096x4096 pixels each, covering a FOV of 5 3 ' . 
The instrument is sensitive to the wavelength range 0.8 - 2.5/jm, 
thus covering the I, Y, J, H, K bands. For astrometric experiments 
the use of the data analysis software Astro-WISE (Valentijn et al. 
(^007 1) is foreseen. 

In order to achieve its science goals (see Sect. 2 for details), 
MICADO needs to reach a stable (time scales of years) astromet- 
ric accuracy of approximately SOyuas. At present 8-lOm class tele- 
scopes, accuracies of a;0.5% of a resolution element can be reached 
regularly (e.g. Fritz et al. (2009)). Therefore from simple scaling 
of results our goal a priori appears reasonable. However, at levels 
of the order of 100//as there are several sources of statistical and 
systematic errors which need to be taken into account carefully. In 
this article we discuss those effects and analyse strategies to bypass 
them. We conclude that reaching an astrometric accuracy of better 
than 50pas is highly challenging in terms of instrument design and 
data calibration but feasible. 

MICADO's astrometric performance should be of the same 
magnitude as that of the future astrometry space mission GAIA 
(e.g. Jordan ((2008)). GAIA will achieve accuracies better than 
a;50yuas only for bright (V<15.5) targets and only at the end of 
its mission. MICADO is expected to achieve this accuracy for tar- 
gets with /^AB < 26. Other space missions like SIM PlanetQuest 
(e.g. Edberg et al. ( [2007] )^ or JASMINE (e.g. Gouda et al. ( f2007t ) 
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also aim specifically at bright targets in order to reach accuracies 
of » 10/jas (at best). 

For illustration purposes, Fig.[T]shows simulated observations 
of the nuclear star cluster of the Milky Way using both present day 
8-lOm class telescopes and E-ELT/MICADO. Physical parameters 
of the star cluster (stellar density profile, luminosity function) are 
taken from Genzel et al. J2003l >. We discuss technical details of our 
simulations in Sect. 4.1. These maps demonstrate the impressive 
progress to be expected with MICADO. 

Although this study is set up for the specific case of MICADO, 
most of its results are valid in general and therefore of interest be- 
yond the E-ELT community. 

This paper is organised as follows. In Section 2, we discuss 
the science cases identified for MICADO. In Section 3, we review 
the concepts and techniques of accurate astrometry. In Section 4, 
we identify and analyse sources of systematic errors one by one 
and describe methods for minimizing those errors. We provide a 
summary of our results and an overall error budget in Section 5 and 
present our conclusions in Section 6. 



2 SCIENCE CASES 

As part of the instrument design study, the MICADO collaboration 
has identified and analysed (Renzini et al. ( 12008 1 . and references 
therein) several science cases for which the high astrometric accu- 
racy of E-ELT/MICADO is crucial and promises major discoveries. 
We discuss them in the following one-by-one. 

2.1 Galactic Centers 

Located at a distance of a;8 kpc, the nuclear region of the Milky 
Way is the closest galactic nucleus, hosting the supermassive black 
hole (M. a 4 X 10^ Mo) Sgr A* (e.g. Gillessen et al. ( [2009] l). It 
is therefore a unique laboratory for exploring the regime of strong 
gravity, accretion onto black holes, and the co-evolution of dense 
star clusters and active galactic nuclei. 

Present-day NIR instrumentation, e.g. VLT/NACO, provides 
astrometric accuracies down to kO.3 mas and angular resolutions 
down to w50 mas (e.g. Fritz et al. ([2009)). This allowed to identify 
several stars on Keplerian orbits around Sgr A* with orbital periods 
down to a;15 years and pericenter distances as small as a 100 AU 
(«12 mas on sky; e.g. Gillessen et al. (2009) ). It made possible to 
study in detail the kinematics and the composition of the nuclear 
star cluster in the gravitational potential of the central black hole. 
With E-ELT/MICADO one can expect to achieve sensitivities that 
are more than five magnitudes fainter than for VLT/NACO. Angu- 
lar resolutions and astrometric accuracies should also improve by 
factors of about five, meaning that proper motions of order 10//as/yr 
(400 m/s) can be detected within few years of observations. Such 
instrumental performance is necessary in order to adress several 
new questions (e.g. Gillessen et al. ( 12009 l l): 

• Identification of stars on closeby Keplerian orbits with periods 
of few years. 

• Measuring the prograde relativistic orbit precession and test- 
ing other effects of general relativity. 

• Probing possible retrograde orbit precession due to an ex- 
tended mass component built from compact stellar remnants. 

• Analyzing the various separate kinematic structures of the nu- 
clear star cluster, and searching for new ones. 

• Quantifying the binary star fraction in the nuclear cluster. 




Figure 1. An illustration of the expected performance of the E- 
ELT/MICADO system. These simulated maps show the central 1" X 1" (i.e. 
8000AUx8000AU) of the nuclear star cluster of the Milky Way at 2.2//m. 
Top panel: The target region as observed with present day 8- 10m class tele- 
scopes. The diffraction-limited resolution is aiSOmas. For comparison with 
actual observations, see, e.g., Genzel et al. (2003), Ghez et al. (2005). Bot- 
tom panel: The same field as seen by MICADO. The angular resolution is 
ailOmas. The improvement in detail and depth is obvious. 



With MICADO, this type of analysis can be extended to other 
nearby galaxies. One obvious example is the core of M31 which 
hosts a M. « 1.4 X 10| black hole. Similar to the case of the Milky 
Way, M3rs nucleus shows several distinct stellar populations: a 
triple nucleus and two nested star disks around the central black 
hole have been identified (Bender et al. ( 12005 1 ). Although M31 is 
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more distant from earth by a factor » 100 compared to the Galactic 
center, the larger mass of its black hole (by a factor of x;35) causes 
stellar proper motions of about 6% compared to those in the nu- 
clear cluster of the Milky Way. Therefore kinematic analyses ana- 
loguous to the Galactic center experiment will require capabilities 
as predicted for E-ELT/MICADO. 

Another example is Centaurus A which hosts a M. a 5 x lOg 
black hole (Neumayer et al. ( 2007 1). Given its distance of a:5 Mpc, 
proper motions of a;10yuas/yr correspond to a;200 km/s. Thus mea- 
suring the motions of stars in the vicinity of the black hole is pos- 
sible. Other galaxies might be interesting targets as well. 

2.2 Intermediate Mass Black Holes 

The expected high astrometric accuracy of MICADO opens a new 
window in the search for and analysis of intermediate mass black 
holes (IMBH), objects with masses of few thousand solar masses. 
In the last years, the detection of those objects in the Arches clus- 
ter (Portegies Zwart et al. d2006) ), the star association GCIRS13 
in the Galactic center (e.g. Maillard et al. ( 2004)), ai Cen (Noyola 
et al. ( I2008t ). and other locations has been claimed. Most of these 
analyses are based on radial velocity dispersion profiles. This in- 
troduces systematic ambiguities as anisotropic velocity dispersions 
can mimick the presence of central point masses; this one can see 
for example in the anisotropy term in the Jeans equation. Therefore 
reliable (non)detections of IMBHs require measurements of stellar 
proper motions (e.g. Anderson & van der Marel ( |2009> ). Typical 
velocity dispersions a, of star clusters are of the order 10 km/s. 
This corresponds to a;50/jas/yr at a distance of 40 kpc, meaning 
that for most of the Galactic star clusters a proper motion analysis 
is feasible with MICADO only. This allows 

• contraining black hole masses in Galactic star clusters within 
few years, 

• probing the low-mass end of the M. - cr, relation, 

• testing the dynamical evolution of star clusters. 

2.3 Globular Clusters 

As discussed above, the astrometric accuracy of MICADO should 
allow measuring stellar proper motions of few km/s for most Galac- 
tic globular clusters. Additionally, direct measurements of cluster 
parallaxes become possible: for a distance of 40 kpc, the full paral- 
lax displacement is 50;uas, corresponding to MICADO's predicted 
astrometric accuracy. This allows adressing several topics: 

• The spatial distribution of globular clusters. 

• Cluster proper motions and their orbits around the Milky Way 
(e.g. Bedin et al. ( [2003112006) ). 

• Internal cluster kinematics, including rotation. 

• Separating cluster members from field stars, thus making anal- 
yses of cluster star populations more reliable (see also Anderson et 
al. ( |2006> ). 

2.4 Dark Matter In Dwarf Spheroidal Galaxies 

Cold dark matter models predict high mass densities and cuspy 
density profiles for the central regions of galaxy halos. In contrast, 
warm dark matter models predict substantially lower central densi- 
ties and constant density cores at small radii. The dwarf spheroidal 
satellite galaxies of the Milky Way provide a unique laboratory to 
test those models. Their proximity of a;100 kpc makes it possible 



to resolve individual stars and analyze their dynamics in the gravi- 
tational potentials of their galaxies. Present day studies are usually 
based on line-of-sight velocity dispersion profiles. However, there 
are degeneracies of velocity dispersion anisotropies with mass den- 
sity profiles. Therefore any conclusive analysis requires measuring 
all three components of the velocity vectors of the tracer stars (e.g. 
Strigari et al. ( l2007t ). 

This type of studies requires accurate proper motion measure- 
ments with uncertainties of few km/s or better. MICADO will be 
able to provide accuracies of «5 km/s within few years of observa- 
tions for targets about 100 kpc away. This makes MICADO a deci- 
sive tool for testing the validity of present-day dark matter models. 



3 THE ASTROMETRY PROBLEM 

Throughout this paper, we use the term "astrometry" the following 
way. We discuss time-resolved relative positions, meaning the po- 
sitions of a science target with respect to a set of reference sources. 
Science target and reference sources are located in the field of view 
(FOV) of the camera, i.e. they are recorded simultaneously in the 
same science image. 

Intra-epoch measurements cover timelines that are so short 
that intrinsic motions of science targets or reference sources cannot 
be detected. This can be a set of images taken within the same night 
or a few adjacent nights. Any variations in measured positions are 
due to measurement errors and can be used to determine the po- 
sition accuracies. All information obtained from this data set (i.e. 
images, coordinates, ...) can be combined, e.g. for improving the 
signal-to-noise ratio (SNR). 

Inter-epoch measurements cover timelines sufficient for de- 
tecting source motions. Those data cannot be combined in a straight 
forward manner. Intrinsic source motions and measurement errors 
interfere and need to be disentangled, usually meaning that cal- 
culating errors requires additional information. Given MICADO's 
proper motion accuracies of a;50yuas/yr or better (depending on the 
duration of the experiments), corresponding to ^10 km/s at 40 kpc 
distance, intrinsic motion of Galactic stars is detected easily. There- 
fore it will be necessary to use extragalactic sources, including high 
redshift objects like QSOs, as references for some science cases. 

If a set of science images is at hand, the general analysis recipe 
is as follows. 

Step 1. From each image «, one extracts the detector posi- 
tions jX") of all sources of interest (science targets and reference 
sources). Detector positions need to be measured with high accura- 
cies of order few milli-pixels (mpix). Existing centroiding, source 
profile fitting, and point spread function (PSF) correlation algo- 
rithms provide such accuracies (e.g. Diolaiti et al. ( I2000t : Berry 
& Bumell (2000 1; Trippe ( 2008 1). 

Step 2. For each image n, the detector positions jX") need to 
be converted into global astrometric coordinates {x" j. The reference 
frame can be the detector coordinates of a selected zero-point im- 
age or any more general astrometric coordinate system. Using the 
detector positions of the reference sources (X"^j) and their astro- 
metric positions {x"^j), one calculates a transformation 

■■ (x;;,) — > {x;;,! . (i) 

Obtaining the positions {x"^jj requires some prior knowledge on the 
reference sources. One powerful approach is cross-calibration with 
other datasets, maybe from other wavelength regimes. A nice ex- 
ample is given by Reid et al. ( (2007T l. They use precise VLBI astrom- 
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etry of SiO maser stars in order to define an astrometric reference 
frame in K-band images of the Galactic center. 

If the reference sources are known to be not moving (e.g. ex- 
tragalactic objects), one can set {x"^j) = {X"); the index indicates 
the selected zero-point in time. The same can be done if the ref- 
erence source ensemble is (or is defined as being) at rest in aver- 
age, i.e. ({xj'j.f)) = <{x"|.f)) (e.g. a sufficiently large set of stars in a 
star cluster). In this case however, one will loose information on 
a global motion (drift, rotation, contraction, etc.) of the combined 
system "science target + reference sources". In any case, the trans- 
formation T„ is used to compute global astrometric coordinates for 
the science targets like 

{X"l — . T„ ({X"l) = {x"l . (2) 

Commonly, low-order (<5) 2-dimensional polynomial coordinate 
transformations 

x' - aQ + aix + a2y + ojX^ + a^y^ + a^xy + ... (3) 
>■' = bQ + bix + b2y + biX^ + Z?4_v^ + bjxy + ... (4) 

are used. In case prior knowledge on the geometry of the required 
transformation is available, one can use models with smaller num- 
bers of free parameters (e.g. Montenbruck & Pfleger ( |1989t ; An- 
derson et al. C2006J ; Trippe et al. (2008 1). 

The number of available reference sources governs the maxi- 
mum order of coordinate transformations (see step 2). A 1st order 
polynomial transform with six parameters requires three reference 
sources, i.e. 2x3 = 6 coordinates. A 2nd-order transform (12 pa- 
rameters) requires six reference sources, and so on. Throughout this 
paper we assume that any two science images need to be connected 
via full astrometric transformations. This means that we regard 
more simple methods of data combination like stacking, simple- 
shift-and-add etc. as non-astrometric and thus not usable for our 
purpose. 

In order to judge the astrometric accuracy achievable with a 
system like MICADO, one has to distinguish statistical and sys- 
tematic influences. The statistical measurement accuracy is given 
by 

A I I A \/5m\/ 100 \ 

Here A is the wavelength, D the telescope aperture, and SNR the 
signal-to-noise ratio (Lindegren (1978)). To give an example for 
the case of the E-ELT: with A = 2.2yum (K-band), D = 42m, 
and SNR = 100, one obtains a statistical astrometric accuracy 
(Ti = 34/jas. From this we see that - in principle - astrometric 
accuracies of ~10/jas can be obtained with the E-ELT. This means 
that any source of additional, especially systematic error needs to 
be compensated down to this level if one actually wants to fully 
exploit the E-ELT's capabilities. 



4 ERROR SOURCES 

In total, we have identified ten effects that might have the potential 
to reduce the expected astrometric accuracy of MICADO substan- 
tially. We will discuss these "Terrible Ten" in the following subsec- 
tions. The first three phenomena we analyze are instrumental, the 
next five are atmospheric, and the last two are astronomical. 



Table 1. Strehl ratio estimates as obtained from preliminary simulations 
(Liske 2008; M. Kissler-Patig priv. comm.). The /iccnier are the central wave- 
lengths of the filters. 



Band 


I 


J 


H 


K 


'iccnlcr [/Jm] 


0.900 


1.215 


1.654 


2.179 


Strehl ratio [%] 


2 


18 


35 


53 



4.1 Sampling and Pixel Scales 

The detector position of a point source can be computed only if the 
source PSF is sufficiently sampled. If the pixel scale - expressed in 
angular units per pixel - is too large (undersampling), position in- 
formation is lost because there is no unique mathematical descrip- 
tion for the PSF profile anymore. Especially, a PSF can then be 
modelled by profiles with different centers of light (i.e. different 
detector positions). This effect is known as the pixel phase error. It 
can reach magnitudes of several tenths of a pixel, thus providing an 
important boundary condition for the instrument design. A detailed 
description of this phenomenon is provided by Anderson & King 
( [2000t . 

In order to identify the critical pixel scale of MICADO, we 
first created artificial PSFs for each of the bands I, J, H, and K. We 
modelled each PSF P(x) as a superposition of a 2-dimensional Airy 
function A(x) and a 2D Moffat profile M(x) like 

P(x) = aA(x) + (1 - a)M(x) . (6) 

For each PSF, a stochastic optimization routin^l adjusted the pa- 
rameter a e [0, 1] such that the resulting PSF profile showed the 
proper Strehl ratio. We took Strehl ratio estimates for MICADO 
from preliminary E-ELT adaptive optics (AO) system simulations 
(Liske (2008) ; M. Kissler-Patig priv. comm. ), the values are shown 
in Table[T] 

For each of the four wavelength bands we examined pixel 
scales from 1 to 7 mas/pix in steps of 0.1 mas/pix. In each con- 
figuration, we placed the PSF at a random detector position and 
re-binned it to the corresponding pixel scale. After this, we fit the 
detector position with a 2D Gaussian light distribution. We iterated 
this procedure 250 times for each configuration. Pixel phase errors 
were the rms values of the distributions of the differences between 
true and measured positions. 

The results of our analysis are shown in Fig. |2] For clarity, 
we restrict the diagram to the decisive pixel scale range from 2 
to 4 mas/pix. For all wavelength bands the errors are smaller than 
sil/jas for pixel scales below 3 mas/pix. Given that the filter bands 
span a factor of 2.4 in wavelength and thus in A/D, it might not 
be obvious why the critical pixel scales we find are that similar 
for all bands. This effect is caused by the differences in Strehl ra- 
tios: whereas the diffraction limited profile width decreases with 
decreasing wavelength, the Strehl ratio decreases, too, meaning 
a stronger atmospheric blur. In our case, the two efi^ects roughly 
counterbalance each other; we actually find the smallest beam size 
for the H band PSF due to its combination of high SNR (a;35%) 
and small diffraction limit (A/D a; 8. 1 mas). This makes it possible 
to quote a single critical pixel scale for all bands. We therefore con- 
clude that for isolated point sources pixel scales up to 3 mas/pix can 

^ Implemented in the MPE data processing soft- 
ware DPUSER developed by Thomas Ott; see 
|http : //www .mpe .mpg . de/~ott/dpuser/history . html| 
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pixel scale [mas/pix] 

Figure 2. Pixel phase eiTor vs. pixel scale for I, J, H, K bands for isolated 
sources. Model PSFs ai'e supeipositions of Airy and Moffat profiles with 
Strehl ratios as given in Table[T] For all bands, the errors are below ailyuas 
for pixel scales smaller than 3 mas/pix. 



be used for MIC ADO without introducing noticable (alpas) pixel 
phase errors. Therefore we use this scale as the standard value for 
the MICADO design. 

In case of crowding, the astrometric accuracy is additionally 
limited by source overlap. For this reason the MICADO design 
foresees a "small scale mode". As the camera design is catoptric 
and all mirrors are fixed, we use a pick-up arm to image a small (6", 
using one out of the 4x4 detectors) FOV with a reduced pixel scale 
of 1.5 mas/pix. The reduced scale helps to better separate close by 
sources and such reduce the misplacements of the source centroids. 

In order to quantify this effect, we examined our two pixel 
scales, 1.5 mas/pix and 3 mas/pix, for each of the bands H and 
K, meaning four configurations in total. For each configuration we 
simulated a crowded star field by placing PSFs at random positions. 
We used the same model PSFs as for the case of isolated sources. 
We applied a luminosity profile corresponding to the K-band lu- 
minosity function of the Galactic bulge which is dlogn/dlogS ~ 
-0.8 (Zoccali et al. J2003I I; n is the source number, 5 is the flux). 
The dynamic range of the star sample was 10 magnitudes, i.e. a fac- 
tor of 10,000 in flux. Source densities were k3000 stars per square 
arcsecond. We also used these routines (with modified parameters) 
to create the map shown in Fig.[T] 

In each simulated map, we searched for stars and calculated 
their detector positions using the PSF fitting routine StarFinder 
(Diolaiti et al. (20001). We derived median uncertainties vs. fluxes 
from the distributions of the differences between true and measured 
positions. As the absolute values of these errors are functions of 
several parameters like fluxes, luminosity profiles, and source den- 
sities, we converted our results into relative numbers, using the er- 
rors obtained for the 3 mas/pix scales as references. 

The outcome of our analysis is presented in Fig. [3] Reducing 
the pixel scales from 3 mas/pix to 1.5 mas/pix improves the typical 
accuracies by factors of a 2-3. The effect is stronger in H than in 
K band; this indicates that images with smaller beam sizes profit 
more from a small pixel scale. For the brightest sources (more than 
a;500 units) the differences between the pixel scales are (at least 
in K band) not very pronounced any more because (a) the number 
of very bright sources is small and (b) they outshine most of their 
neighbours. We can thus conclude that for crowded fields the use 
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Figure 3. Relative median astrometric error vs. source flux for crowded 
sources, separately for H and K bands. The results for the pixel scales 
1 .5 mas/pix are given as fractions of the uncertainties found for 3 mas/pix. 
The horizontal dashed line indicates a ratio of 1. Error bars indicate the 
68% uncertainty ranges. This analysis shows that reducing the pixel scale 
improves the typical accuracies by factors sj2-3. 



of small pixel scales down to 1.5 mas/pix is indeed important for 
keeping a high level of astrometric accuracy. 

Combining all results for isolated and crowded sources, we 
can conclude the following for the design and the operation of MI- 
CADO: 

• For sufficiently isolated target sources, pixel scales smaller 
than or equal 3 mas/pix are free of noticable («lyuas) systematic 
uncertainties and suited for accurate astrometry. This scale is there- 
fore going to be the standard pixel scale for MICADO. 

• For the special case of crowded sources, smaller pixel scales 
down to 1 mas/pix can substantially (factors a 2) reduce the as- 
trometric errors introduced by source overlap. We therefore im- 
plemented a "small-scale mode" with a reduced pixel scale of 
1.5 mas/pix into the MICADO design. This mode will be used to 
map very crowded regions. 



4.2 Instrumental Distortion 

The geometric distortion of an optical system can seriously limit 
its astrometric accuracy. In the following, we discuss non-linear 
distortions. Linear terms like shifts, rotations, scalings, or shear, 
are absorbed by 1st (or higher) order coordinate transforms (see 
Sect. 3). In case of non-linear distortion, the effective pixel scale 
is a function of detector position (e.g. Greason et al. ( il994i l. and 
references therein). 

For MICADO, the amount of distortion to be expected (mean- 
ing the difference between imaged and theoretical image positions) 
is of the order of few per cent. The largest numbers, about 5%, we 
find for the case of an Offner design for the camera. Although other 
optical designs might have the potential to reduce the amount of 
distortion down to «0.3% (Dierickx ( |2008» , it cannot be neglected. 
Across a FOV of 1', a distortion of 0.3% corresponds to a misplace- 
ment of 0.18". For comparison: an astrometric accuracy of a;50yuas 
across the same FOV corresponds to a relative error of less than 
10"*. Such low uncertainties cannot be achieved by design. In or- 
der to reach the desired astrometric accuracies, the distortion must 
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be compensated by dedicated calibration schemes - regardless of 
its amplitude. 

Another effect to be considered is imperfect fabrication of the 
detectors. The HAWAII-4RG detectors to be used for MICADO are 
designed with pixels sizes of 15;um. The position accuracy of these 
pixels in the detector grid is a lOnm; for a pixel scale of 3 mas/pix 
this corresponds to an astrometric uncertainty of a; 2/jas (Richard 
BlanlH priv. comm), i.e. is hardly noticable. However, any distor- 
tion calibration scheme must provide the ability to catch potential 
inaccuracies of the detectors (see also Anderson ( 12002) for the case 
of HST/WFC). 

Depending on the complexity of the distortion, two methods 
for its description are possible. Analytic descriptions make use of 
analytic parametrizations analogous to the coordinate transforms 
discussed in Sect. 3. Most commonly, 2D polynomials up to about 
5th order are used as distortion models (e.g. Greason et al. ( |1994| l, 
and references therein). This approach allows covering all effects 
up to a selected order without prior knowledge on the geometry of 
the problem. The disadvantage of this method is the large number 
of model parameters to be calculated. A recent example for this 
approach is the analysis of the nuclear star cluster of the Milky Way 
by Ghez et al. (20081 and Lu et al. ( 20091. They use a polynomial 
model to correct the geometric distortion of the NIRC and NIRC2 
cameras at the W.M. Keck Observatory. 

If a physical model for the distortion is at hand, a correspond- 
ing model can be substantially simpler. A recent example for this 
approach is the analysis of the nuclear star cluster of the Milky Way 
by Trippe et al. ( |20U8] l and Gillessen et al. (( |2009] l). They use the 
3rd-order model 

r = r'(l-/?r'') (7) 
with 

r = X - Xc and r' = x' - 

(e.g. Jahne ( |2005t fl in order to correct the distortion of the imager 
NAOS/CONICA at the VLT. Here x and x' are the true and dis- 
torted image coordinates respectively, y6 is a parameter describing 
the strength of the grid curvature, and x^ = {xc,yc) is the zero 
point of the distortion on the detector. This approach has the obvi- 
ous advantage that it requires only three parameters {j3,Xc,yc) for 
a 3rd-order distortion model compared to 20 for the case of a full 
polynomial solution. The most important disadvantage is the need 
for accurate a priori knowledge of the distortion geometry. 

In cases where analytic solutions are not feasible or not ac- 
curate enough, empiric descriptions might be used instead or in 
addition. This means that the information of interest is stored in 
look-up tables. Those tables usually have the dimensions of the de- 
tector(s) and give the amount of distortion (or correction) for each 
pixel (separately for x,y). This approach is necessary if significant 
high-frequency distortion - meaning that the spatial scales of the 
signal are small compared to the detector size - is present. For 
the specific case of MICADO, the gaps between the detectors will 
introduce discontinuities into any astrometric solution. This effect 
could be caught by using lookup tables. 

In a detailed analysis of the WFC instrument aboard the Hub- 
ble Space Telescope, Anderson ( 12002) uses a combined "polyno- 

' Teledyne Imaging Sensors, Camaiillo, USA 

See also the electronic manual of the public Gemini North Galactic Cen- 
ter Demonstration Science Data Set for another application on Galactic 
Center imaging data. 



mial model plus lookup-table" ansatz to model the distortion of the 
camera. He achieves residuals smaller than 0.01 pixels with this 
method. For the case of MICADO, this would mean residual errors 
below «30yuas which is an amount acceptable for our purpose. 

There are several methods to extract the distortion parameters. 
On-sky methods make use of dedicated observations of astronom- 
ical targets; usually, star clusters are used (e.g., Anderson ( 2002)). 
If the true (astrometric) coordinates of the sources in the target field 
are known, one can derive the distortion parameters by comparing 
true and observed (detector) positions. If the true source positions 
are not known, one can observe the target field many times with 
slightly different telescope pointings. In this case, one compares 
the pairwise distances of objects that are present in two or more 
images. Modulations in these distances as function of detector po- 
sition are equivalent to (non-linear) distortion. 

By construction, on-sky methods are sensitive to distortions 
introduced by the atmosphere. Given the high accuracies we seek, 
those methods are not practicable (at least stand-alone) for the case 
of MICADO but will be a secondary approach. 

In-lab methods characterize the instrument with dedicated 
measurements in the laboratory or at the telescope. One example 
is the "north-south test" used for the spectro-imager SINFONI at 
the VLT (e.g. Abuter et al. (2006)). These methods use devices that 
illuminate the detectors with well-defined images or light patterns. 
Comparing the theoretical with the observed images allows for a 
description of the distortion. 

Given that in our acuracy regime the atmosphere can severely 
limit the quality of on-sky calibration images (see Sect. 4.4 - 4.8 
for details), we assume that we need to implement an internal cal- 
ibration device into MICADO. For this purpose, we examined the 
use of a calibration mask located in the focal plane of the imager. 
Such a mask could be a regular pattern of holes in an intransparent 
material. However, the following calculation shows that this is ac- 
tually challenging. With MICADO, we want to observe a FOV with 
extension x = V with an accuracy 6x = \Qp.as, meaning a relative 
accuracy of 

6x1 XK 1.7 X 10"^ . 

The width of MICADO's focal plane will be approximately / = 
0.25m. This scales to a positioning accuracy for any reference lo- 
cated in the focal plane of 

6l = lx Sxjx » 4.2 X 10"'*m . 

As we see here, the positions of the holes in our calibration mask 
need be known with accuracies of about 40nm (the hole diame- 
ters would be a:30/im corresponding to diffraction-limited point 
sources). This result does not mean that the fabrication process 
needs to be accurate at this level. Instead, it is sufficient to map 
the hole positions with accuracies of s;40nm after the making of 
the mask. State-of-the-art photolithographic techniques provide the 
required production and mapping accuracies at standard tempera- 
tures (a:300K). With this approach, the calibration mask would be 
a transparent quartz plate with a chrome cover into which the holes 
are etched. Although such a system is thermally very stable, the im- 
pact of possible non-linear thermal deformation needs to be inves- 
tigated via dedicated laboratory experiments (B. Lorenfl D. Ros^l 
priv. comm.). We consider the calibration mask approach to be the 

^ Center for NanoSciences, Ludwig Maximilians University, Munich, Ger- 
many 

* Rose Fotomasken, Bergisch Gladbach. Germany 
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easiest one in terms of telescope and instrument design. We iiave 
not yet investigated in detail more "exotic" ideas, e.g. illuminating 
the detectors with a well-defined diffraction pattern or spectrum. 

As we see from the discussion in this section, it is crucial to 
gather the maximum amount of information on the instrumental 
distortion of MICADO. We therefore conclude that the following 
steps are necessary: 

• Estimating amount and geometrical structure of the distortion 
theoretically from the optics design. 

• Careful mapping of the camera in the lab and at the telescope. 
Implementation of an internal calibration device, e.g. a calibration 
mask. 

• Testing the system on-telescope for any evolution of the dis- 
tortion. Evolution parameters can be time (aging effects), telescope 
orientation (gravitational flexure), etc. 

• Additional regular dedicated on-sky calibration observations 
of sufficient astronomical targets, e.g. star clusters, as secondary 
tests. 



4.3 Telescope Instabilities 

Instabilities of the telescope system have the potential to affect as- 
trometric experiments. For the present design of the E-ELT one ex- 
pects relative intra-night plate scale variations of a;0.1% (Gonzalez 
( 12008 ^). Across a FOV of 1', this corresponds to position variations 
of order 60 mas. Another effect are adapter-rotator instabilities that 
can introduce systematic frame-to-frame rotation. Those rotations 
introduce position misplacements of the same order of magnitude 
(see, e.g., Trippe et al. ( 12008) for the case of VLT/NACO). 

Fortunately, those effects - shifts, global scalings, rotations - 
are linear in geometry. Therefore they can be controlled via coordi- 
nate transforms of 1st or higher order without additional calibration 
steps. This statement does not hold however for gravitational flex- 
ure effects that introduce a time-variable non-linear instrumental 
distortion. This phenomenon is covered by our discussion of in- 
strumental distortion in the previous subsection. 



4.4 Achromatic Diflferential Atmospheric Refraction 

Any ground-based position measurement is affected by atmo- 
spheric refraction. As we discuss relative position measurements 
of sources located in our FOV, we have to take into account any 
relative or differential atmospheric refraction. In this case, we have 
to discriminate achromatic and chromatic effects. 

Achromatic differential refraction is caused by the slight dif- 
ference in zenith angles ( of two (or more) sources located within 
the same FOV. For each source, the atmospheric refraction leads to 
a deviation between physical and observed zenith angles. For two 
sources at slightly different zenith angles, those deviations will be 
different. Therefore the observed distance between the two targets 
deviates from the physical one and needs to be corrected. The re- 
quired correction Ax of the distance between two sources and 1 
is approximately given by the relation 



Ax = (l+ tan^ (i)(A + SStan" ^i)Af . 



(8) 



Here is the zenith angle of source 1, A<f is the observed zenith 
separation, A and B are constants. Detailed calculations (Gubler 
& Tytler ( |1998l l) show that the linear terms of this effect are of 
the order of several milli-arcsec, whereas the quadratic terms are 
as small as ^lyuas. For example: with = 45° and Af = 30", the 
correction amounts to » 15 mas in the first-order terms and a;2//as in 



second order. Fortunately, these terms can be absorbed by quadratic 
coordinate transforms (see Sect. 3). From this we conclude that we 
can neglect achromatic differential atmospheric refraction for the 
purpose of our analysis. 



4.5 Chromatic DilTerential Atmospheric Refraction 

Chromatic differential refraction (CDR) or atmospheric dispersion 
is a more severe problem for accurate astrometry than the achro- 
matic case. As the refractive index n of the atmosphere is a function 
of wavelength, the observed angular distance between two sources 
is a function of the (relative) source colours. For a given true zenith 
distance the deviation from the apparent zenith distance ^„ (in 
radians) follows the approximative relation 

<r, - f„ - Rtan(, = ('Lzl\tan(, 



2n 



(9) 



Here R is the refraction constant. For standard conditions the re- 
fractive index is given by 



(n-r)x 10" = 64.328 -l- 



29498.1 X 10- 



255.4 X lO^" 

+ .. , T (10) 



146 X 10-« - i- 41 X 10-'^ - s 



with s = /i"', A being the vacuum wavelength in nm (e.g. Cox 
( I2000t ). Of course, this has impact on relative astrometry only if 
the two sources have different colours. For observations of two stars 
with broadband JHK filters, astrometric errors are of the order of 
1 mas (within a wide range, depending on zenith angles, angular 
distances, and source colours). We therefore need to correct the 
CDR in order to meet our desired accuracies. 

In the following, we investigate the use of an atmospheric dis- 
persion corrector (ADC) placed into the optical path. This ADC is a 
pair of ZnSe/ZnS biprisms which refracts - for a given zenith angle 
- the infalling radiation such that the CDR is compensated. In order 
to use this element for a range of zenith angles, the biprisms can be 
rotated around the optical axis relative to each other; this controls 
the strength of refraction. Such a system was built and operated 
successfully for the SHARP II-l- NIR camera system (Eisenhauer 
(1998)). The required optical components are available in the spa- 
tial dimensions needed for MICADO (S. Koebele priv. comm.). In 
order to decide on the design strategy, we explored three options in 
detail: 

(i) The use of one non-tuneable ADC optimized for simultane- 
ous correction of the full wavelength band range J, H, K, (A = 
1.1...2.35yum). 

(ii) One ADC covering the entire band range that can be tuned 
in wavelength to bands J, H, and K,, respectively. 

(iii) Three ADCs, one for each of the three bands J, H, and K,. 

For each option (1,2, 3), we computed residual dispersions (disper- 
sions left after correction) vs. wavelength for each of the bands J, 
H, and K,. All calculations used zenith angles ^ = 45°. Our compu- 
tations made use of numerical optimization routines implemented 
into the software package Mathematic^ The results for J band are 
given in Fig.|4] Whereas option 1 results in a very asymmetric curve 
of correction, options 2 and 3 are very similar in shape and highly 
symmetric. As we see here, one can reduce the residuals by using 
an ADC specifically designed for the wavelength band analyzed 
(option 3) by factors of about ten compared to a tunable ADC (op- 
tion 2). However, the impact on astrometry is actually given by the 
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Table 2. Relative zenith angle shifts computed for the options, science 
cases, and filter bands discussed in the text. All numbers are in /jas. 



Band 


J 


H 


Ks 


Option 1: 








Case A 


730 


73 


196 


CaseB 


642 


112 


250 


Option 2: 








Case A 


9 


9 


6 


CaseB 


35 


23 


21 


Option 3: 








Case A 


9 


9 


2 


Case B 


35 


22 


11 




1.2 



1.3 



X [/.im] 



Figure 4. Residual dispersion vs. wavelength in J band after correction with 
an ADC assuming ( = 45°. Top panel: Using one non-tunable ADC for the 
entire JHK band range. Center panel: Using one tunable ADC for the entire 
JHK band range tuned to J band. Bottom panel: Using an ADC specifi- 
cally designed for optimum correction in J band. Please note the changes 
in scales. The non-tunable ADC results in a very asymmetric correction 
curve. The other two curves are offset by about ten times their amplitudes 
from each other, but are very symmetiic in their shapes. 



relative displacement of two sources with different colours. This 
means that the residual astrometric error is mainly controlled by 
the symmetry of the curves, not by their amplitudes or absolute 
levels. We may therefore expect that options 2 and 3 result in very 
similar astrometric accuracies. 

For a quantitative description of the residual astrometric errors 
we analyzed two hypothetical science cases. We examined 

• an observation of two black bodies of very different tempera- 



tures (r, =3,000 K, Tn =30,000 K), crudely corresponding to M5 
and BO main sequence stars, respectively ("case A"), 

• the case of two black bodies at T =5,800 K (e.g. sun-like stars) 
affected by very different extinctions of Ay = 25 and Ay = 35, 
respectively ("case B"). For this case, we used the extinction law 
by Draine (1989 ) as taken from Lutz et al. ( 199 6 1. 

For each option, science case, and filter band we computed 
relative shifts in zenith angle between the sources 1 and 2 via 



(11) 



/ S,idA j S,idA 

We present the results of our analysis in Table |2] As we see here, 
using a non-tuneable ADC (option 1) for the entire JHK band range 
leaves us with errors of roughly a lOOyuas in H band which contains 
the center of correction. When going to Kj and J bands, the resid- 
uals increase to a;200/jas and s;700/ias, respectively. All numbers 
exceed any acceptable value by factors of about ten or more. This 
clearly rules out option 1 . 

Using a tuneable ADC (option 2) or three ADCs optimized 
for the three bands (option 3) reduces the residual zenith angle dis- 
placements to a;10/ias. This result complies with our requirements 
and shows that our approach can indeed achieve the necessary ac- 
curacies for realistic science cases. As expected from Eqs.|9land|10l 
the errors increase with increasing frequency. Therefore one should 
avoid observations at very short (shorter than J) wavelengths. Ad- 
ditionally, one should consider the use of narrow-band filters in 
case of (a) observations at short wavelengths or (b) extreme rela- 
tive source colours. 

Another calibration step one should consider is an a posteriori 
correction. As we discuss above, the impact of CDR on astromet- 
ric solutions can be quantified analytically if the relative source 
colours are known. However, this requires careful monitoring of 
the atmosphere (see, e.g., Helminiak (2009)). 

From our analysis of chromatic differential atmospheric re- 
fraction, we conclude the following: 

• CDR distorts astrometric solutions by up to few mas depend- 
ing on source colours. We need to implement a dedicated correc- 
tion. 

• Using a tuneable ZnS/ZnSe atmospheric dispersion corrector 
reduces the astrometric signal caused by CDR to a;10/ias in JHK 
bands asuming realistic science cases. We therefore highly recom- 
mend to implement such a device into MICADO. 

• In case of extreme relative source colours or observations at 
short (shorter than J-band) wavelengths, one should consider the 
use of narrowband filters additional to an ADC. 



High Precision Astrometry with E-ELT/MICADO 9 



• If the source colours and the atmospheric conditions (temper- 
atures, pressures, humidities) are known with sufficient accuracies, 
one might additionally apply an analytic a posteriori correction of 
astrometric data. 



4.6 Guide Star Measurement Errors 

MICADO's MCAO system is designed to make use of three natu- 
ral guide stars that are located (at a priori arbitrary positions) in the 
field of view. Additionally, the use of six laser guide stars is fore- 
seen. The stars are observed simultaneously as references for cor- 
recting the wavefront deformation imposed on astronomical signals 
by the atmosphere. The natural guide stars are used for correcting 
low-order effects, whereas the laser guide stars correct high-order 
distortions. The LGS high-order correction does not make use of 
the guide star positons. 

This is different for the NGS low-order correction. This cali- 
bration step requires knowledge of the relative guide star positions. 
The guide star measurements will have finite errors due to atmo- 
spheric fluctuations that introduce a "position wander". Uncertain- 
ties in the measured guide star position will introduce time- variable 
distortion into the AO corrected FOV. For N stars, one may expect 
distortions up to order N-1, i.e. up to 2nd order for the case of 
three natural guide stars. These distortions will differ from image 
to image. Therefore it is necessary to combine images with full co- 
ordinate transform of minimum order two. If this is done, the effect 
will be compensated completely. 



4.7 Differential Tilt Jitter 

The light from a science target and the light from an adaptive op- 
tics system reference source travel through different columns of at- 
mospheric turbulence. An AO system applies a tip-tilt correction to 
the signal received from the AO reference source. This correction is 
slightly different for other positions in the field of view. Therefore, 
any two objects in the observed field suffer from differential tilt jit- 
ter: a random, achromatic, anisotropic fluctuation of the observed 
angular distance of the two sources (e.g., Britton J2006I I; Cameron 
et al. J2009l >). In first order, the rms of this fluctuation follows the 
relation 



ctt J oc ex D 



-7/6 



(12) 



Here 9 is the angular distance between the sources, D is the telecope 
aperture, r is the aperture wind crossing time (approximately: D 
divided by wind speed), and / is the integration time. Tilt jitter is an 
anisotropic effect with 

V3o-_L « 1.732o-_L 



(13) 



where iT||, cr± denote the tilt jitter rms parallel and perpendicular to 
the line connecting the two sources, respectively. 

In Fig. [5] we give an example for a tilt jitter signal observed 
in images taken with an 8-m-class telescope. This result is from 
the work by Fritz J2009I I who analyzed diffraction-limited VLT im- 
ages of the nuclear cluster of the Milky Way. For pairs of stars, he 
computed the uncertainties of the measured distances parallel and 
perpendicular to the lines conecting the two stars. Fig.|5]shows the 
uncertainty as function of star angular distance. By means of linear 
fits one finds a relation crn/(T± = 1.91 ± 0.22, i.e. V3 within errors 
- as expected fom Eq.[T3l From the line slopes and the known in- 
tegration times one can estimate an aperture wind crossing time of 
T X 0.6s. For details refer to Fritz l [2009]t and Fritz et al. ( |2009l l. 



O parallel 
■ perpendicular 




["] 

Figure 5. A differential tilt jitter signal as observed in VLT images. This 
diagram shows the en'or on the measured distance between two stars o" as 
a function of star angular distance &. Errors are calculated parallel (||) and 
perpendicular (±) to the line connecting the two stars on sky. Points with 
errorbars are data, continuous lines indicate the best fitting linear models. 
From the line slopes one finds o"||/(T^ = 1.91 ± 0.22, i.e. V3 within errors 
as expected fom Eg. 1131 This result is taken from Fritz (2009). 



The impact of this effect on astrometry can be substantial. To 
give a reference: based on observations with the Hale 200 inch 
telescope, Cameron et al. J2009i ) find astrometric uncertainties of 
a;75 mas for D = 5m, 9 = 60", and / = t = 0.2s (i.e. r/f = 1). 
For a realistic infra-epoch E-ELT observation with D = 42m and 
t = IOOt k 100s (i.e. r/t = 0.01), this scales to cttj ~600y[(as. Av- 
eraging out this error to reasonable scales (less than 50/jas) would 
take about four hours. This imposes severe constraints on any as- 
trometric observation. 

The calculation given above is strictly valid only for single- 
conjugate AO (SCAO) systems, i.e. AO systems using one ref- 
erence source. However, MICADO is designed to use an MCAO 
system with three natural guide stars initially. The use of multiple 
guide star should reduce the tilt jitter error substantially. Ellerbroek 
( 2007 1 finds that for a 30m telescope and the case of three guide 
stars arranged in an equilateral triangle the error is reduced by a 
factor a;6 compared to the SCAO case. Scaling the results by Eller- 
broek {2Mf\ to the case of the E-ELT, one finds 

(Ttj « 430 X r"^ 



■ pas 



(14) 



with t being the integration time in seconds and 9 = 60". This 
means that accuracies of a;10/ias can be achieved within integration 
times of about 30 minutes. 

A helpful property of tilt jitter is the fact that it is random in 
time, but correlated in space. This is also reported by Fritz l |2009l > 
who finds that already using linear coordinate transforms between 
images reduces the error by a factor a;2. According to Cameron 
et al. (2009), the use of coordinate reference frames based on 
weighted pairwise distances between a target source and several 
reference sources is able to catch the tilt jitter error. This approach 
is similar to the concept of kriging (e.g. Clark & Harper ( 2000)). 
It requires some tens of astrometric reference points in the FOV in 
order to obtain a sufficient number of pairwise baselines. Addition- 
ally, it requires any intra-epoch dataset to consist of some tens of 
individual exposures instead of few long-term integrations. This is 
necessary in order to estimate the uncertainty (and thus the weight) 
of each pairwise source distance from its histogram. If this scheme 
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Figure 6. Difference image of two simulated MAORY K-band PSFs located 
70" away from each other at opposite sides of the center of AO correction. 
This image covers a field of » 7.5" X 6". The difference map shows a 
complicated, highly symmetric pattern. 



is applied, the error can be reduced to «10/jas within few minutes 
of integration. 

From our discussion on differential tilt jitter we conclude for 
MICADO: 

• Differential tilt jitter can introduce errors of order alOOyuas 
into typical E-ELT/MICADO observations. This effect is not a chal- 
lenge for the instrument design but for the planning of observations 
and data calibration. 

• The tilt jitter error integrates out with time as oc r''^. For the 
E-ELT, reducing the error to order of a;10/jas, integration times of 
at least a;30 min per intra-epoch dataset will be necessary. 

• The tilt jitter error can be calibrated out (down to a;10ytias) by 
using dedicated astrometric reference frames and transforms. This 
requires (a) some tens of reference sources in the FOV and (b) some 
tens of individual exposures per intra-epoch dataset. 



4.8 Anisoplanatism 

All existing AO systems suffer from anisoplanatism, meaning that 
the shape of a PSF is a function of its position in the field of view. 
Usually, this spatial variability is described by the relation 



(5r [(Lias] 

Figure 7. Astrometric errors introduced by MCAO anisoplanatism in H {top 
panel) and K (bottom panel) bands. Errors are centroid shifts (with respect 
to the case of perfectly symmetric light distributions) of the difference im- 
ages between pairs of simulated PSFs located at different locations in the 
MAORY FOV. The histograms have cutoffs at xlpaf, in H band and at aSpas 
in K band. 



impact of anisoplanatism is a systematic degradation of the PSF's 
Strehl ratio S like 



5(0) = 5oexp[-(e/0„>^/3] 



(16) 



P{g) = p(0) * K{0) 



(15) 



(Fusco et al. J2000l l; Steinbring et al. J2002l |2005t ; Cresci et al. 
( 12005 ^). Here is the angular distance from the center of AO cor- 
rection, P{Q) [P(dy\ is the PSF at the center of correction [at dis- 
tance 0\, and K{8) is a kernel describing the PSF variation; * de- 
notes convolution. 

If K{6) were asymmetric, then the center of light of the PSF 
- and thus the detector position - would be a function of 6. This 
would introduce a systematic distortion of astrometric solutions 
that has to be calibrated out. 

For existing SCAO systems, K{6) is known to be symmetric. 
It can be approximated analytically as an elliptical Gaussian pro- 
file (e.g. Steinbring et al. l l2005b ). In this case, the most important 



assuming a perfect AO correction. Here is the Strehl ratio at the 
center of AO correction, the parameter 6n is the anisopanatic an- 
gle. This effect reduces the SNR of a source and thus the statistical 
astrometric accuracy, but a priori does not introduce a systematic 
error. 

For the case of MCAO systems with multiple centers of cor- 
rection the problem is less well understood. In order to estimate 
the astrometric error introduced by anisoplanatism, we made use 
of preliminary simulations of multi-guide-star MAORY PSFs in H 
and K bando The simulation data provide a grid of 217 PSFs lo- 
cated at distances between 0" and 80" from the center of correc- 
tion. The simulations assume a seeing of 0.6". PSFs are sampled 
with pixel scales equal to A/2D, corresponding to 4.05 mas/pix and 
5.30 mas/pix for H and K bands, respectively. 



^ Made available by the MAORY 
|http : //www . bo . astro . it/~maory/Maory/| 



consortium 
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By construction, each PSF is centered accurately at the cen- 
tral pixel of a 512x512 pixel grid. As discussed above, we need to 
estimate the relative position error introduced by anisoplanatism. 
For this we computed for each pair of PSFs a difference map by 
subtracting one PSF image from the other In Fig. |6] we show the 
difference image for the case of two PSFs located 70" away from 
each other at positions (-60", 0") and (+10", 0") with the center of 
AO correction located at (0, 0). Clearly, the map shows a compli- 
cated, highly symmetric light distribution. From the high degree of 
symmetry one can already suspect qualitatively that the impact of 
MCAO anisoplanatism on the astrometry is small. 

In order to quantify this statement, we calculated for each dif- 
ference map the centroid 

2,- ■ Ij 

^centroid — ~. (17) 

of the light distribution. Here x is an arbitrary coordinate, x, and 
/, are position and flux assigned to the i-th pixel. The summation 
is performed over all map pixels (e.g. Berry & Bumell (2000^). As 
we had 217 PSF images at hand, this procedure resulted in 23436 
measurements. Centroid positions different from the center of the 
pixel grid correspond to position shifts of one PSF with respect to 
the other caused by asymmetries in the anisoplanatism kernel. 

We summarize our results in Fig.|7] The histograms of the po- 
sition shifts show quite sharp cutoffs at a;7yuas in H band and at 
~8//as in K band. We do not find a significant correlation of posi- 
tion errors with distance from the center of correction or any other 
systematic relation. 

We conclude that MCAO anisoplanatism probably introduces 
very small but noticable uncertainties into astrometry. Depending 
on wavelength band we find errors (upper ends of histograms) up 
to aSyuas. As the anisoplanatism effect is not (or very weakly) cor- 
related with PSF positions in the FOV, it can probably not be caught 
by coordinate transforms or any other systematic correction but has 
to be included into the error budget. 



4.9 High-z Galaxies As Reference Sources 

Throughout this article we discuss relative astrometry, meaning 
measuring relative source positions. So far we implicitly assumed 
the reference sources to be point sources, especially stars. How- 
ever, in some cases the use of reference points other than stars will 
be necessary: 

• The number of stars located in the FOV can be too small. As 
discussed in Sect. 3, the number of reference sources is a func- 
tion of the order of the coordinate transform used. For example, in 
case of a 2nd-order transform, the minimum number is six. Using 
a larger number than the theoretical minimum is useful to average 
out statistical position errors that could propagate into the trans- 
formation. Additionally, systematic effects like instrumental distor- 
tion (Sect. 4.2) and differential tilt jitter (Sect. 4.7) require large 
numbers of reference points for calibration purposes. Deep fields 
at moderate galactic latitudes {\b\ » 50°) include only a few to a 
few tens A'ab <23 stars/arcmin (see e.g. Forster Schreiber et al. 
d2006t : Wuyts et al. (2008)). This number is probably too small for 
astrometric MICADO observations. 

• The science case might require an extragalactic reference 
frame, thus excluding Galactic stars. For example, several groups 
(e.g. Bedin et al. ( |2003ll2006l l: Kalirai et al. ( |20D4l l) use samples of 
background galaxies as reference in order to measure the absolute 



motion of globular clusters. Other science cases may need similar 
approaches (see Sect. 2). 

In order to investigate the use of galaxies as references, we 
analyzed two simulated MICADO deep field K-band images as- 
trometrically. We simulated galaxies following realistic distribu- 
tions of K-band magnitude as function of redshift z taken from the 
Chandra Deep Field South FIREWORKS ("field A"; Wuyts et al. 
f2008< ) and FIRES Hubble Deep Field South ("field B"; Labbe et 
al. ( I2003t ) catalogs. All objects in the fields are built using four 
main ingredients: 

• Smooth elliptical and disk galaxy light distributions with real- 
istic surface number densities, sizes, and K-z distributions. The K-z 
distributions were obtained by drawing at random galaxies from 
the two source catalogues above, scaling the numbers by the ratio 
of the MICADO FOV and the respective survey areas. A fixed mix 
of ellipticals and disks was adopted (40%/60%), and effective radii 
were assigned at random so as to match approximately the size dis- 
tributions of Franx et al. ( 12008b . Uniform distributions are assumed 
for the axis ratios (with minimum of 0.3 for ellipticals and 0.1 for 
disks) and for the position angles. 

• Bulges (for disk galaxies) with de Vaucouleurs luminosity 
profiles adopted for simplicity. Bulge properties are poorly con- 
strained at high redshift; our assumptions were guided by results 
from Elmegreen et al. ((2009 1, and references therein; see also Gen- 
zel et al. (i2008i )). We assumed uniform distribution of axis ratios 
between 0.7 and 1, position angles fixed at those of the host disks, 
elfective radii inversely proportional to 1 + z and such that = 2kpc 
at z = 0, and a random distribution of bulge-to-total light ratios be- 
tween and 0.5. 

• Clumps (for disk galaxies) with realistic numbers per disk, 
typical light fractions, characteristic sizes and dependence on red- 
shift, and radial distribution across the disks so as to roughly 
match observed properties at z » 1-2.5 (e.g., Genzel et al. (2008J; 
Elmegreen et al. ( 120091 ) and references therein) and their lack at 
z « 0. Of order 1-10 clumps per disk were simulated, with Gaus- 
sian light profiles (again for simplicity), a narrow distribution of 
FWHMs varying as (1 + z) such that FWHM = 1 kpc at z = 2.5, 
light fractions of a few percent typically, an exponential distribu- 
tion of radial positions with scale-length five times that of the host 
disk, random position angles and axis ratios between 0.7 and 1. 

• Unresolved star clusters (for all galaxies). This is the ingre- 
dient which is least constrained at high redshift, as current instru- 
ments lack the resolution and point-source sensitivity to detect indi- 
vidual clusters if such exist at high z. To define their properties, we 
assumed superpositions of five dilTerent populations of point-like 
sources with a range of (Gaussian) magnitude and radial distribu- 
tions. The brightest population contains the fewest clusters, with 
the most centrally concentrated distributions. The brightest cluster 
in each galaxy contributes 0.1% of the total light. These assump- 
tions roughly reproduce those of super star clusters and globular 
clusters in local galaxies, and unresolved sources in local analogs 
of z 3 Lyman-break galaxies (Overzier et al. ( I2008t ). 

For each field image we created three realizations corresponding to 
E-EUT/MICADO integration times of 1, 4, and 10 hours, respec- 
tively, by adding Gaussian noise to the original map. As given by 
the input catalogs, field A has a point source brightness limit of 
^AB = 24.3, field B of ATab = 25.6. These limits are 5it depths 
for point sources in circular apertures of 2" diameter. In Fig. [8] we 
show four simulated galaxies from field B as examples. 
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Figure 8. A few simulated galaxies from our artificial deep field B. Images 
in the top line cover x 3.5" X 3.5", images in the bottom line cover » 
1.8" X 1.8". 



The star clusters are especially interesting for astrometry. 
They provide an ensemble of several thousands of bright, point- 
like sources that can be analyzed with standard tools developed for 
point source astrometry. In each field we selected all star clusters 
that were sufficiently isolated, meaning a distance of at least 30 mas 
(about three resolution elements) to the next source. In total we 
used 1600 sources in field A and 2174 sources in field B. In each 
image we fit the clusters with 2D Gaussian brightness distributions 
in order to accurately (meaning few milli-pixels in the best cases) 
determine their positions. We then compared the results with the 
true source positions on the original noise-free maps. The distribu- 
tions of the deviations between true and measured positions provide 
measures of astrometric accuracies. 

We show the results for lOh integration times in Fig. [9l For 
both fields we find distributions with Icr- widths of a;920yuas (field 
A) and s;650/ias (field B) per coordinate. These numbers corre- 
spond to the typical measurement accuracies for individual star 
clusters. The global astrometric accuracy is given by the standard 
error of the mean of the distribution, i.e. cr/ Viv with A' being the 
number of sources. 

We summarize the astrometric accuracies found for all fields 
and integration times in Table [3] Not surprising, we find the high- 
est accuracies of 23yuas and 14/jas (per coordinate) for fields A and 
B, respectively, for the largest integration time of lOh. The errors 
increase to about 65/jas and 29ytias at integration times of Ih. The 
systematic difi'erence between fields A and B originates from their 
diff'erent depths: field B is based on a catalog that has a » 1 .3 magni- 
tudes deeper brightness limit; it is therefore "richer" of sufficiently 
bright targets. As MICADO will perform astrometric observations 
of sources A"ab < 26 (see also Sect. 1), the results found from 
field B are actually more realistic for describing the performance of 
MICADO. We can therefore conclude that total (i.e. the quadratic 
sum of the values for both coordinates) astrometric accuracies of 
x:20;uas are realistic for integration times of order lOh when using 
galactic star clusters as referencess. 
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Figure 9. Histograms of position errors found from using several thousand 
point-like star clusters as references. These results are for fields A (top 
panel) and B (bottom panel) at lOh integration time. Please note the dif- 
ferent number axes scales. 



So far, we have taken into account non-resolved (point-like) 
sources as references only. However, the global extended light 
distribution contains additional information. We therefore com- 
puted 2D cross-correlations between different realizations of the 
same "observations". As reference we used one 7000x7000 pixels 
(21"x2r') sub-field of field A containing seven galaxies. For each 
integration time (Ih, 4h, lOh) we computed three realizations of the 
random noise map added to the original light distribution. For all 
pairwise combinations we computed the 2D cross-correlations. We 
calculated their centers by means of 2D Gaussian fits to the cen- 
tral parts of the cross-correlation maps. The uncertainty of the map 
centers provides a measure of astrometric accuracy. For all integra- 
tion times we found very similar results: 4.5/jas at Ih, 4.1/jas at 4h, 
and 3.9yuas at lOh. This indicates that at these level the accuracy is 
only weakly correlated with the SNR but dominated by systematic 
effects. In any case the impressive accuracies of a:4/jas suggest that 
using 2D cross-correlations of images can provide highly accurate 
astrometry. 

However, 2D cross-correlations cannot be used for relative 
astrometry in a straight forward manner. By construction, cross- 
correlations are only sensitive to shifts between images. They can- 
not describe transformations like scalings, rotations, or higher order 
transforms as discussed in Sect. 3. Using them for astrometry there- 
fore requires a more sophisticated approach than computing the 
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Table 3. Astrometric accuracies found from using several thousand point- 
like star clusters as references. We give errors for the coordinates x and y as 
function of field and integration time. 
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cross-correlation between two (or more) images. One ansatz might 
be to calculate cross-correlations for a sufficient number (about ten 
or more) of sub-fields of the target area. This would provide a set of 
relative position (changes). Feeding this dataset into a proper model 
might result in a highly accurate coordinate transform. 

From our analysis we conclude that the use of high-z galax- 
ies as astrometric calibrators is feasible at the level of a;20yuas at 
integration times of about lOh. This number assumes the use of 
several thousand point-like star clusters as reference sources (also 
assuming that such exist at intermediate to high redshift). Cross- 
correlation of images might provide even better accuracies, but 
would require a sophisticated computational scheme that can ac- 
tually provide full coordinate transforms. 

In some cases, one might consider separating the intra-epoch 
and the inter-epoch calibration steps. Firstly, one can use a small 
number of stars (as discussed earlier, six may be sufficient) for in- 
tra-epoch (frame-to-frame) calibration. Proper adding of images al- 
lows building up SNR sufficient for detection of and position mea- 
surements on faint galaxies. Secondly, one can use these galaxies 
for inter-epoch calibration. 



As discussed in Section 2, the science cases for E- 
ELT/MICADO demand accurate measurements of positions and 
proper motions over spatial scales from few ten yuas (e.g. parallaxes 
of globular clusters; see Sect. 2.3) to few hundred mas (e.g. stellar 
orbits around Sgr A*; see Sect. 2. 1). A relative pixel scale accuracy 
of 2 X 10"^ corresponds to an error of 10 yuas over an angular dis- 
tance of 500 mas. From this we take the message that the accuracies 
of the reference positions should not exceed a 1 mas substantially. 

There are several possibilities for obtaining very accurate ref- 
erence source positions. One of them is the use of sources located in 
the MICADO FOV and visible in both NIR and radio, e.g. QSOs or 
maser stars. If radio-interferometric positions (e.g. VLBI) - which 
are accurate on the sub-mas level - are at hand, the uncertainty on 
the pixel scale can be very low (see, e.g., Reid et al. ( I2007t for 
the case of the Galactic center). Another option might be the use 
of stars from the HIPPARCOS catalogue (Perryman et al. ( 11997) ). 
However, the uncertainties in stellar proper motions, which are 
s;0.8 mas/yr, will propagate into position uncertainties of a:20 mas 
over a timeline of 25 years (given that the HIPPARCOS catalogue 
reference epoch is J1991.25; Perryman et al. ( |1997| l). This means 
that the astrometric uncertainties will exceed 10 yuas if the spatial 
scale of the experiment is larger than a; 25 mas. A more promising 
approach is the use of data from the future GAIA astrometry space 
mission that will provide absolute astrometric accuracies well be- 
low 1 mas (Jordan ( |2008l l). The main advantage of the GAIA cat- 
alogue (its high accuracies aside) is the large number of stars in- 
cluded - «10' - meaning that for any arbitrary field in the sky a 
sufficient number of reference points should be available. 

We conclude from our discussion that for MICADO sufficient 
numbers of calibration sources with absolute position accuracies 
better than 1 mas will be available. We therefore expect that the 
calibration of the projected pixel scale introduces errors of ~10/ias 
at most. 



4.10 Calibration of the Projected Pixel Scale 

In Section 3 we have outlined the concepts of relative astrometry 
and coordinate transforms. These schemes make use of sets of ref- 
erence positions {x"^^) in order to compute transformation matrices 
for the «-th dataset. However, all transformation and calibration 
steps we have discussed up to now are executed in image space, 
i.e. in units of pixels. After removing non-linear distortions from 
the data, one needs to calibrate the pixel scale as projected on sky 
in order to accurately convert measured positions and motions into 
angular units. We note that the projected pixel scales for the two 
coordinates x, y may be dilferent if the detector plane is tilted with 
respect to the focal plane. 

Calculating the scaling factors requires astrometric reference 
points with known positions located in the target FOV. The number 
of reference sources should be at least three in order to allow for 
a full linear transformation. Fortunately, this calibration step needs 
to be executed only once for a selected "master" (or "zero") image. 
All other images can be connected to the master image reference 
frame via coordinate transforms, including the proper scaling (e.g. 
Trippe et al. (2008)). Errors on the reference positions propagate 
into the positions and motions calculated from the data. For the 
FOV of MICADO, reference position accuracies of order 1 mas 
translate into a relative scaling accuracy of 

6x1 lmas/53" « 2x 10"^ 



5 RESULTS AND ERROR BUDGET 

We have identified and discussed ten elfects that might limit the ex- 
pected astrometric accuracy of E-ELT/MICADO observations sys- 
tematically. We have been able to quantify each of these sources of 
error. From this, we can calculate a prediction for the error budget 
of MICADO. 

(i) For isolated point sources, detector sampling / binning does 
not introduce noticable (a;lyuas) errors as long as the pixel scale 
does not exceed 3 mas/pix. For sources affected by crowding, using 
a smaller scale of 1.5 mas/pix can improve the accuracy by factors 
of about two compared to the 3 mas/pix case. Therefore MICADO 
will use a pixel scale of 3 mas/pix as standard and a reduced scale 
of 1.5 mas/pix for mapping crowded fields. For the error budget, 
we can thus assume a sampling error 

0"samp = 1 yUaS . 

(ii) Instrumental geometric distortion needs to be taken into ac- 
count by dedicated calibration procedures. We propose to imple- 
ment a calibration mask into MICADO that illuminates the de- 
tectors with a well-defined image. Such a mask would have to be 
mapped with accuracies of «40nm. Based on our results in combi- 
nation with published works using Hubble Space Telescope data, 
we estimate that distortion can be corrected down to levels of » 10- 
30/ias. For the error budget, we therefore use 
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o-jis, = 30 //as . 

(iii) Telescope instabilities, notably plate scale instabilities and 
instrumental rotations, are linear effects that can be absorbed by 
coordinate transforms. Therefore they do not contribute to the error 
budget. 

(iv) Atmospheric achromatic differential refraction is important 
(order 10 mas) only in linear terms which can absorbed by coordi- 
nate tranforms. Higher-order contributions are of order lyjas, mean- 
ing for the error budget 

o"ADR = 1 /ias . 

(v) Chromatic differential refraction introduces position errors 
of the order al mas in NIR observations depending on relative 
source colours. A tuneable ZnS/ZnSe atmospheric dispersion cor- 
rector can reduce this effect to al0~20yuas for "typical" science 
cases. Extreme relative source colours might require the additional 
use of narrowband filters and/or analytic a posteriori correction 
schemes. For the error budget, we thus set 

o-cDR = 20 yuas . 

(vi) AO guide star measurement errors for N natural guide stars 
can introduce distortions up to order N- 1 into images, meaning 2nd 
order distortions for the three NGS of MICADO. This effect can be 
absorbed by coordinate transforms (of order N- 1). It therefore does 
not contribute to the error budget. 

(vii) Atmospheric differential tilt jitter can introduce errors of 
alOO/jas into diffraction-limited E-ELT observations. It integrates 
out with r"-. For MICADO which uses an MCAO system, the tilt 
jitter error can be integrated down to «10/jas within about 30 min- 
utes of observation. Using dedicated coordinate transform allows 
reaching this accuracy in shorter times. For the error budget, we 
thus set 

cttj = 10 yuas . 

(viii) The anisoplanatism of the MAORY AO system introduces 
uncertainties of up to a;8/ias. There appears to be no correlation of 
this error with the position of a PSF in the FOV; it therefore can not 
be calibrated out in a straight forward manner. Therefore we add it 
to the error budget: 

o"ai,iso = 8 /jas . 

(ix) Depending on target field and science case, the use of galax- 
ies as astrometric calibrators may be necessary. From a simulated 
MICADO deep field we find that we can use several thousand non- 
resolved galactic star clusters as point-like reference sources. How- 
ever, good accuracies of a:20/jas require long intra-epoch integra- 
tion times of about 10 hours; we consider this to be a somewhat 
large but realistic timescale. We therefore use for the error budget 

0"galaxics = 20 /jaS . 

(x) The accuracy of the sky-projected pixel scale is limited by 
the accuracy of astrometric positions of reference sources in the 
FOV. Given the typical accuracies of present-day radiointerfero- 



metric data and the future GAIA catalogue, which are better than 
=s 1 mas, we add to the error budget 

fT-scale = 10 ^aS . 



From the individual uncertainties listed above we can calcu- 
late a total intrinsic astrometric accuracy for MICADO as 




This number provides a systematic limit for astrometric accuracies 
to be expected from MICADO data. Of course, this result corre- 
sponds to a somewhat arbitrary "typical" case. As many parameters 
like integration times, source colours, numbers and types of refer- 
ence sources, etc. can vary over wide ranges, the actual iXsys for 
a specific observation can be quite different - in both directions - 
from the one we quote here. Nevertheless, we conclude that we are 
able to quantify the mean systematic astrometric accuracy achiev- 
able with MICADO which is a;40yuas. 

When discussing the accuracy of the measurement for a given 
target, one of course needs to add the statistical measurement error 
o-L (Eq.O which scales with the SNR. For SNR=100, ctl = 34/jas, 

and thus the combined error is .^tr^y, -I- cr^ = 56pas. SNRs differ- 
ent from 100 modify this number accordingly. 



6 CONCLUSIONS 

In this article we have studied the capabilities expected for the NIR 
imager MICADO for the future 42-m European Extremely Large 
Telescope with respect to accurate astrometry. A variety of science 
cases requires long-term astrometric accuracies of «50/jas. We dis- 
cuss and quantify ten effects that potentially limit the astrometric 
accuracy of MICADO. We conclude that the systematic accuracy 
limit for astrometric observations with MICADO is o-^y^ « 40yuas. 
We find that astrometry at this accuracy level with MICADO re- 
quires the fulfillment of several conditions: 

• All images, regardless of their distance in time, need to be 
combined via full coordinate transforms of second order or higher. 

• MICADO needs to be equipped with an astrometric calibra- 
tion mask for monitoring the instrumental distortion. The pixel 
scale of the camera should not exceed the 3 mas/pix used in the 
current design. 

• Astrometric observations require decent integration times of at 
least 30 minutes per epoch. This is unavoidable in order to average 
out atmospheric tilt jitter. When using high-z galaxies as astromet- 
ric reference points, integration times up to about 10 hours can be 
necessary. 

It is noteworthy that the effects discussed in this article already 
affect observations collected with present 8m-class telescopes. In 
his exhaustive analysis of NIR images obtained with VLT/NACO, 
Fritz ( 2009 J has been able to detect signatures of chromatic differ- 
ential refraction and differential tilt jitter in his astrometric dataset. 
He concludes that taking into account these effects can improve 
the accuracies down to few hundred yuas. This agrees with the find- 
ings of Lazorenko ( 12006b and Lazorenko et al. ( 12007 l l who analyze 
seeing-limited optical R-band (/Iccntcr = 655nm) images taken with 
VLT/FORSl+2. They conclude that they are able to achieve astro- 
metric precisions (but not accuracies) of a lOOyuas by using a special 
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scheme for scheduling observations and dedicated coordinate trans- 
forms (although they neglect instrumental geometric distortion). 

The analysis we provide here is set up for the specific case 
of E-ELT/MICADO, but parts of our results are valid in general. 
This study should thus contain valuable information for other fu- 
ture 30-40m telescopes. As some of the effects we discuss are ac- 
tually observed in present day 8m-class telescope data, our analysis 
might also be helpful for the calibration of data already taken. We 
therefore expect that our work is of interest well beyond the E-ELT 
community. 
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